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ABSTRACT 

It is shown that above the temperature of maximum abundance, recombination rates 
into the excited states of Hc-hke ions that are calculated using earlier, more approx- 
imate methods differ markedly from rates obtained from recent distorted-wave and 
R-Matrix calculations (unified recombination rate coefficients) for Ca, Fe, and Ni. 
The present rates lead to G ratios that are greatly lower than those resulting from the 
more approximate rates in previous works, by up to a factor of six at high electron tem- 
peratures. Excellent agreement between the distortcd-wave and the R-Matrix rates, as 
well as excellent agreement in the G ratios calculated from them, provides support for 
the accuracy of these new values which have a broad applicability to the modelling and 
interpreting of X-ray spectra from a variety of astrophysical and laboratory sources. 

Key words: atomic data - atomic processes - X-Rays: general ~ line: formation 



1 INTRODUCTION 

The Ka emission lines of helium-like ions ha ve been actively 
studie d for quite some time, starting with lEdlen fc T^oinl 
Much of the basic theory and interpretation that 
is now used for astrophysica l spectroscopy was established 
by lGabriel fc JordanI (|l969al ) in their work to interpret so- 
lar spectra. While detailed discussions of the structure of 



Gabriel & JordanI Il969al 


bl iGabriell 


I972I: iBlumenthal et al.l 


I972I: Mewe & Schriivei 


Il978a bl: 


Pradhan & ShuUl 198ll: 


Pradhan 19821"). in brief there are four He-like lines in the 



Is" 'So 



Ka complex: ls2p 'Pi Is" 'So (w), ls2p ^P2 
(x), ls2p ^Pi ^ Is" 'So (y), ls2s =*Si Is" 'So (z). The 'w' 
line is sometimes referred to as the resonance or dipole al- 
lowed line and 'y' as the intercombination line as it is dipole 
allowed through fine-structure mixing. The lines 'x' and 'z' 
both arise from higher order transitions, and as such are 
described by much slower rates. The 'x' line represents a 
magnetic quadrupole transiti on, and the 'z' lin e repres ents 
a magnetic dipole transition. iGabriel fc JordanI l|l969ah de- 
fined the temperature sensitive line ratio G as 



G = 



J(x) + J(y) + /(z) 
J(w) 



(1) 



where /(w) is the total emission in the 'w' line in units of 
number of photons cm~^ s~^. In general, the behaviour of 
the G ratio is determined by the rate at which the He-like 
excited states are populated. At high temperatures these 
populations are determined by an interplay between the pro- 
cesses of recombination and impact excitation. Recombina- 
tion preferentially feeds the three lines arising from triplet 
states ('x', 'y', and 'z') while electron-impact excitation pref- 
erentially feeds the 'w' line. If the recombination rate is dom- 
inant, due, for example, to an increase in the number of H- 
like ions relative to the number of He-like ions, one would 
expect to see more flux in the 'x', 'y' and 'z' lines than in the 
' w' line (as is observed in the case of recom bining plasmas 
(|Pradhanlll985l : IOelgoetz fc PradhanI 120041 ) ). resulting in a 
G ratio greater than one. 

For heavier atoms the situation is more com- 
plex as satellite lines produced by radiative decay 
from autoionizing three-electron states are a prominent 
part of the Ka spectrum (see Gabriel fc JordanI 1969bl: 
Gabriellll972l: iMewe fc Schriiveilll978al: 1 Swartz fc SulkanenI 
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I993I : iBautista fc Kallmanll2000l : Toelgoetz fc PradhanlbOOlL 
among others). These satellite lines become less important 
as temperature increases. As the high-temperature regime 
(above 10^ K) is the interest of this present work, no in- 
vestigation of line ratios involving satellite lines has been 
made, although we note that satellite lines can be, and 
often are, important at temperatures well above 10^ K 
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l|Oelgoetz fc PradhanI I2OOII ). If the spectra are not ade- 
quately resolved, these satellite lines should be included in 
any line ratios used for analysis. 

As this study investigates the effect of recombination 
rates (both radiative (RR) and dielectronic (DR)) on the 
high-temperature behaviour of the G ratio, a brief re- 
view of the data commonly used, how they are used, and 
the reported high-temperat ure dependence of the G ra- 
tio is provided. CHIANTI (iDere et alj Il997l : iLandi et all 
I2OO6I). XSTAR jBautista fc Kallmanll200lf l. and ATOMDB 
( Smith et al.l 120011 ) are three commonly cited databases in 
astrophysics. All three use different types of recombination 

data. 

CHIANTI uses the recombination data of I Me we et al.l 

(Il985h for most helium-like ions which is based on the work 
of lBurgess fc Seatoij (Il960l) on recombination into neutral 
he lium. The work of iMewe fc Schri ivcr (T978a) is also based 
on iBurgess fc SeatonI (| 19601 7 and has been used in some re- 
cent investigations such a s lOelgoetz fc PradhanI (I2OOII ) as 
well as the earlier work bv iPradhan fc Shulll ( 198ll ) . Fe and 



Ca are the exceptions to the above rule; for thes e elemen ts 
CHIANTI uses the work by iBelv-Dubau et al.l (Il982allbh 
which is based on earlier studies of h ydrogen by Burgess! 
1 19581 ) . For completeness we note that iMewe et al. ( 19851 ) 
also depart from the above met hodology and inste a d direc tly 
fit recombination data from iBelv-Dubau et al] (|l982al la) 
for Fe XXV and Ca XIX . The hydrogenic methodology of 
iBelv-Dubau et al.l (|l982al lbh has also been used in calcula- 
tions for other elements (e.g. IPorguet fc Dubaul (|2000l )). 

XSTAR takes a different approach; it calculates recom- 
bination rates from the photoio nization cross sections pro- 
duce d by the Opacity Project (jSeatonl Il987l ; I Cunt o et al.l 
1 19931 ). using more approximate methods only for n > 10. 
These Opacity Project cross sections were ca lculated us- 
ing R-Matrix method s in LS couplin g (iBurke ct al. I9Y1I: 
iBerrington et al]|l987l ). Most recentlv lOelgoetz"fc PradhanI 
(|2004 ) used unified recombination rates based on R-Matrix 
calculations (|Nahar fc PradhanI 1 1994 ) with relativistic fine 
structure included in the Breit-Pauli approximation. 

ATOM DB uses the f i ts to distorted-wave photoioniza- 
tion data of IClark et al.l (|l986l ) for He-like cross sections. 
These cross sections are integrated to obtain the corre- 
sponding rates in the standard manner. Additionally there 
are numerous individual calculations that use recombina- 
tion rates from distorted-wave data in a similar fashion. 
Thus, much of the recombination data falls into three ba- 
sic categories: data derived from m ore a pproximate meth- 



P 

ods (s uch as those bv lBurge"s3 (|l958l ) and lBurgess fc Seaton 
( 1960I )). data based on R-M atrix calculations (|Burke et al.l 
I I971I ; iBerrington et al.lll987l ). and data based on distorted- 
wave calcu lations. 

Both iMewe et all (| 19851 ) and jBautista fc KallmanI 

(|2000|)_ calculate the emission intensities of Fe XXV as 
does jPradhan fc Shulll (|l98ll ), using the iMewe fc Schriiveij 
(|l978al ) rate coefficients. If one calculates G ratios from the 
data presented in table IV of iMewe et al.l (|l985l ) the value 
of the G ratio is steadily decreasing, with values of approxi- 
matel y 0.5, in between 10* and 10^ Kelvin. iPradhan fc Shulll 
l|l98ll ) show a G ratio taking a sharp upturn in the same 
temperature range and h aving a value greater than 1.0. 
iBautista fc KallmanI (|2000l ) for this same temperature range 
report a G value in the range of 0.25 for Fe XXV, with a 



possibly slight rise as the electron temperature approaches 
10® K. Thus the disagreement among these calculations is 
significant and resolving this discrepancy i s the g oal of the 
present work. To this end we use GSM (jOelgocta .20061 ) ■ 
a code recently developed to model plasma emission spec- 
tra, to explore the effects of different types of recombina- 
tion data on the G ratio at high temperatures. Data from 
recent distorted-wave and R-Matrix calculations are used 
and compared to rates c alcul ated using expressi ons from 
iMewe fc Schriived (|l978al ) and iMewe et all (|l985l ); the ex- 
act rates in the XSTAR and ATOMDB databases are not 
considered here. Additionally, only the effects of using dif- 
ferent recombination data sets are explored while all other 
rates and the ionisation balance are identical for all models 
of a particular element. 



2 THEORY 

GSM is a code based on the quasi-static approximation, that 
is the excited states of a given ion are considered to always 
be in equilibrium with the ground state of the ionisation 
stage from which they arise, and the neighbouring ionisation 
stages (see lLoch et al.|[2004h . GSM makes the additional ap- 
proximation that some states can be treated purely as con- 
duits that, through radiative decay or auto- ionisation, al- 
low population to transition into the final states of interest, 
thus neglecting the effect of collisions with these intermedi- 
ate states. Both approximations are valid in the case of low 
density plasmas where the radiation field can be neglected 



enced in this work (see 


Gabriel fc Jordan 1969a b; Gabriell 


I972I; Blumenthal et al. 


Il972; Mewe fc Schriivcr 1978a.br 


Pradhan fc ShuU 198ll; Oelgoetz fc Pradhan 2001.. .2004,). 



As a consequence of these approximations the first step 
in any calculation is to solve for the total population of each 
ionisation stage. This is done by solving the following set of 
coupled differential equations: 

dt 

-X,iVe(a,^i-i(e) +C,^,+ i(e)) 



= 7Ve(X,+ia,+i_.,(e)-KX,_iC,-i^,(e)) 



(2) 



where Xi is the total population in the i* ^ ionisation stage, 
A^e the electron number density, e a variable that describes 
the shape of the electron distribution (in the case of this 
work e is the electron temperature as all results presented 
are for thermal systems with a Maxwellian distribution), 
C is a bulk coUisional ionisation rate coefficient, {3 a bulk 
3-body recombination rate coefficient, and a a bulk recom- 
bination rate coefficient (which includes radiative and di- 
electronic recombination). Because all results presented in 
this work neglect the effect of a radiation field, photoioniza- 
tion and stimulated recombination have been omitted from 
equation ([2]). While GSM is capable of calculating these 
bulk rates by integrating the raw level-specific cross sec- 
tions and summing over all intermediate pathways, in gen- 
eral that approach would require a data set that is signifi- 
cantly larger than desired. As such these rates, when avail- 
able, come from either tabula ted or fitted data in the litera - 
ture (such as can be f ound in lArnaud fc Rothenfiu ijl985l ): 
lArnaud fc Ravmondl (|l992h or lMazzotta et all (|l998l )). AU 
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of the calculations presented in this study are steady-state 
calculations and therefore the left-hand side of equation ((2|) 
can be set to zero. We note that GSM is also explicitly de- 
signed to handle transient sources. 

Once the total ionisation stage populations have been 
calculated, they are used to calculate the individual level 
populations, and subsequently the spectra, by solving the 
following set of coupled equations: 



Al-cff 



(3) 



where the variables are defined more or less as before, with 
Ni^j being the population in the j*'' state of the ionisa- 
tion stage, g"*^ an electron-impact (de-)excitation effective 
rate coefficient, an autoionization rate coefficient, 

and A"^^ an effective radiative decay rate. Just as in equa- 
tion (121), equation © omits terms that involve the radiation 
field, ft also omits excitation by proton and a particle im- 
pact, as they are neglected in this study for the sake of sim- 
plicity in order to better illustrate the effect of level-specific 
recombination rates on the G ratio. Once again all calcula- 
tions consider steady-state conditions, so the left-hand side 
of equation ^ can be set to zero. 

As mentioned earlier, GSM does not treat every spec- 
troscopic level the same. Instead, it divides them into two 
categories: explicit (which are included as levels in equation 
(Is}, and coupled via all processes considered), and statisti- 
cal (the effects of which are included when calculating the 
rate coefficients that go into equation Q). The inclusion of 
statistical states in the rate coefficients is accomplished via 
the coUisionless transition matrix (CTM). A given element 
of the CTM can be thought of as the probability that once 
a statistical state gets populated, neglecting collisions, the 
atom or ion will end up in a given explicit state without 
passing through any other explicit states. If Q is the set of 
explicit states, the CTM can be defined using the following 
recursive expression: 



+ 



(4) 



where Ti^k is the appropriate type of spontaneous rate to 
connect i and k, via either radiative decay or autoioniza- 
tion. It should be noted that in order for this quantity to be 
meaningful, the state i must not be in the set Q. 



After calculating the CTM, each effective rate coeffi- 
cient can be constructed from the direct rate coefficient con- 
necting the two explicit states and sums over all the indirect 
paths. For example effective coUisional ionisation rate coef- 
ficients are calculated as 



E 



CfZr{e)Tk- 



(Efc>Ej) 
direct / 



recombination (RR+DR) rate coefficients as 



(5) 



off 



q;c(£) 



RR / \ , r-iDC / \ , 



+ 



E 



D 



E 

^m>E 

DC 



(£)T™_. , (6) 



{Em>Ej ,Em>Ei,E,„>0) 



where is a dielectronic capture rate, and electron- 

impact excitation and de-excitation rate coefficients are cal- 
culated as 



off / N 

gj-*fe(e) 



direct / \ 

Qj^k (e) 



E 



direct/ \m 



+ 



E 



(7) 



(Ei>Ej,Ej>Et.,Ej>0) 



It should be noted that in equations ©-(ITl) the effect of 
resonances are included as indirect pathways via terms in- 
volving dielectronic capture {D^'^). When R-Matrix data 
are used, these terms are omitted to avoid double counting 
the effect of resonances. 

Once equation Q has been solved and the level popu- 
lations obtained, synthetic spectra are calculated. The first 
step is to calculate the intensity (in number of photons 
cm~'^s~^) of all lines in the desired energy range as 



l,k) = Ni^jAj^k 



(8) 



Each line is then given a line shape corresponding to a ther- 
mal Doppler-broadened Gaussian profile. The total spectra, 
S, for a given photon energy, hv, can be expressed as 



S{hv) 



Ishv— e 

2-KkT 



c^(hix-AEij)^ 
2Ae'^ .kT^ 



(9) 



where a ranges over the set of all included transitions in the 
desired energy range, and the ion temperature, Ti, is taken 
to be equal to the electron temperature. 



3 COMPUTATIONS 

The baseline model considered in the present wo r k uses 
the ionisation balance data of iMazzotta et al.l (|l998l ) 
and detailed, level-specific data calculated with the Los 
Alamos suite of atomic p hysics codes (see for instance 
lAbdallah et al.lll994 [ioOlf l. The first step in such a calcu- 
lation is to use the CATS code to calculate the bound and 
autoionizing wave functions, energies, and dipole allowed ra- 
diative decay rates for all fine-structure levels arising from 
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the configurations nl, Isnl, 21nl', Is^nl, ls21nl', and IsSlnl' 
with n ^ 10 and \ ^ g which span the H-like, He-like, and 
Li-like ionisation stages. Only fine-structure levels arising 
from the Is, Is^, ls21, ls^21, and ls21nl' configurations with 
n ^ 10 and 1 ^5 g were treated explicitly. All other lev- 
els were treated statistically. Convergence with respect to n 
was tested with a larger model that included up to n 20. 
The GIPPER code is then used to calculate autoionization 
and photoionization cross sections in the distorted-wave ap- 
proximation, as well as coUisional ionisation cross sections. 
As the ions under consideration are highly charged, a scaled 
hydrogenic approximation was used for coUisional ionisation 
cross sections as it agrees well with distorted- wave results for 
systems such as these and is much more computationally ef- 
ficient. The ACE code was used to calculate distorted- wave 
cross sections for all electron-impact excitation transitions 
out of the lowest seven levels of the helium-like ionisation 
stage, as well as the ls^21 complex of the Li- like ionisation 
stage (which is used to obtain effective coUisional ionisation 
rate coefUcients into the He-like stage as well as a method 
for producing the Ka satellite lines). Lastly, the non-dipole 
A values that give rise to the 'x' and 'z' lines as well as 
a two-photo n decay rate from ls2j_So — > Is'^ ^So were ob- 
tained from iMewe fc Schriiverl (|l978al ). All rates were cal- 
culated from these cross sections inside of GSM. We refer to 
this baseline model as DW. 

The second model, which is designated RM, begins with 
the first model and replaces the radiative and dielectronic 
recombination data from the H-like ionisation stage into the 
He-like ionisation stage with total and level-specific rate co- 
efficients for recombination into all fine-structure levels SLJ 
with n ^ 10 calculate d from Breit-Pauli R-Ma trix photoion- 
ization cross sections (|Nahar fc Pradhan|[l99i ). While there 
are unified recombinatio n data av ailable in the literature for 
Ni (Nahar 200^ 1 and Fe (|Nahar et a l. 20011). we do not con- 
sider a Ca, RM model as there are no data yet available 
for this species. In addition, the energy levels and radiative 
decay rates inside the He-like ion are replaced by R-Matrix 
data (the radiative decay rates ar e taken from Nahar fc Eiss - 
ner (in preparation) for Ni and iNahar fc PradhanI (| 19991 ) 
for Fe). Lastly, the unified recombination rat es through a 
given satellite line of iNahar fc PradhanI (|2006D are used to 
calculate the dielectronic recombination contribution to the 
intensity of all the KLL satellite lines in this model. 

The third model, referred to as M85, also begins 
with the baseline model and replaces the recombination 
rates from the H-like ionisation stage into the He-like ex- 
cited states with rates calculated using the expressions of 
iMewe et all (|l985l l. These expressions are fits to previously 
pu blished data. In the ca s e of N i, the fits are to the rates 
of iMewe fc Gronenschildl (|l98lj ). which are based on the 
work of Burgess fc SeatonI (|l960l ) concerning the recombi- 



nation of neutral He. These fits do not include the con- 
tribution of dielectronic recombination, so in the Ni M85 
model dielectronic recombin ation is included us ing the data 
fr om the baseline D^y model . [Mewe et al.l (Il985l') fit the data 
oflBely-Dubau et al.l (|l982al ') for Fe and iBelv-Dubau et all 
l|l982bl ) for Ca. Both of these works ar e based on the hydro- 
genic reco mbination work of iBurgessI (Il958l). It should be 
noted t hat lBelv-Dubau et al. I (|l982al ) and lBelv -Dubau et al.l 
l|l982bl ) only present rates up t o an electron tempe rature of 
10** K. As such, using the fits of lMewe et"aLl (jigSSl ) to these 



data beyond this temperature may exceed the range of va- 
lidity intended bv lBelv-Dubau et al.l (|l982al lbl). Furthermore 
these rates do include the effect of dielectronic recombina- 
tion, so this contribution is not recalculated and included a 
second time by GSM. 

The fourth model, MS78, is much like the third, but in- 
stead calcul ates radiative recombinati on rates using the ex- 
pressions of lMewe fc Schriivej lll978ah. which are also based 
on the work of lBurgess fc SeatonI I 1960h . As these rates do 
not include dielectronic recombination, this model includes 
DR rates obtained from the data contained in the baseline, 
DW model. 



4 RESULTS 

The results for the G ratio for all models under considera- 
tion are presented in figure[T] The differences are striking. As 
temperature increases and recombination becomes more im- 
p ortant to spectral formatio n, the more approximate rates 
of lMewe fc Schriiverl ([l9783) produce G ratios that diverge 
significantly from the results obtained with the newer DW 
and RM models for all three elements. The M85 result for 
Ni displays similar discrepan cies while the M85 results for 
Fe and Ca (which are fits to iBelv-Dubau et al.l (|l982a| [bh) 
agree better with the RM and DW models. Not surprisingly 
(given the highly charged nature of these systems) the RM 
and DW models are in excellent agreement at all tempera- 
tures. 

Upon investigation of the underlying spectra (figures[2]- 
(4]) to find the source of the divergence in predicted G values 
at high temperature, one continues to see excellent agree- 
ment at all temperatures between the RM and DW models; 
but the analysis is not as straightforward for the more ap- 
proximate methods. In the case of Ni (figure [2} , both the 
M85 and MS78 model agree well with each other, predict- 
ing significantly more fiux in the 'x', 'y', and 'z' lines, and 
slightly more flux in the 'w' line, at high temperatures than 
either of the more accurate models. This behaviour is con- 
sistent with the fac t that both models are ba sed on the same 
underlying work bv lBurgess fc SeatonI (|l960l ). The spectrum 
for Fe (figure[31l is more interesting because, while the MS78 
model continues to predict significantly stronger 'x', 'y', and 
'z' lines, and a slightly stronger 'w' line, compared to the 
DW and RM models, the M85 model predicts a rather dif- 
ferent result. Like the MS78 model, the M85 model predicts 
a slightly stronger 'w' line, but the 'z' line is in much bet- 
ter agreement with the DW and RM models than the MS78 
model. Furthermore the M85 model does a good job of pre- 
dicting the 'x' and 'y' line intensities. 

The detailed results for Ca (figurejj} are somewhat sim- 
ilar to those presented for Fe. The MS78 model continues to 
predict much stronger 'x', 'y', and 'z' lines when compared 
to the DW data set and only a slightly stronger 'w' line. The 
M85 model is in good agreement for the 'x' and 'y' lines and 
much closer for the 'z' line, but in the case of Ca it pre- 
dicts a much stronger 'w' line. This over prediction tends 
to mitigate the effect of the enhanced 'z' line, producing a 
fortuitously lower value for the G ratio. This type of cancel- 
lation underscores the importance of examining the actual 
spectra to gain a better understanding of the quality of the 
underlying atomic data. 
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Figure 1. G ratios for Ni (top), Fc (middle) and Ca (bottom) for each of the models considered. 



In order to determine the fundamental reasons behind 
the above discrepancies in the spectra and G ratios, the re- 
combination rates must be examined in detail, but first it is 
useful to ensure that the models are large enough (as a func- 
tion of the principal quantum number n) to have converged 
upon the correct answer. As the various sources cited in this 
work all chose different limits for the maximum n value used 
in their calculations, lack of convergence could be a reason 
for the differences. Additionally, such a study can be used 
to inv estigate the relative impo rtance of cascade corrections, 
which iPradhan fc ShuUl (|l98ll ) suggested to be the cause of 
the upturn in the G ratio displayed in that work. As the re- 
sults for each element are very similar, only recombination 
rate coefficients for Fe are presented (figure [5} and discussed 
for only the DW model. 



In figure [5] it can be seen that by n = 10 the effective 
rate coefhcients into the relevant four spectroscopic levels 
have converged to within ~5 per cent when the temper- 
ature is above 10* K. All models considered in this work 
extend at least this far in principal quantum number; thus, 
the extent of the model (in n) is not the cause of the discrep- 
ancies. Of additional note is the fact that figure [5] indicates 
that recombination cascades become less important (as a 
method of populating the He-like excited states) at higher 
temperatures, even as recombination is becoming a more im- 
portant process for the population of He-like excited states. 
According to the analysis of the DW model, it is the di- 
rect recombination rate that is the most important process 
at higher temperatures while cascades play a progressively 
smaller role. 



6 J. Oelgoetz et al. 



> 



10.0 
8.0 
6.0 



mi 1 1 III 1 1 II 1 1 II I III 1 1 III 1 1 II I III I III 1 1 III I III III 1 1 III 1 1 III I III III 1 1 III 1 1 III Mill II 1 1 III 1 1 III I III 1 11 1 1 III 1 1 III I III iiM 



a 
> 

^ 4.0 



> 



2.0 

0.0 
0.25 
0.20 



I 0.15 
0.10 



I 0.05 
0.00 & 



> 



0.05 
0.04 



I 0.03 
I 0.02 



0.01 



T =1.0x10 K 



RM 

M85 

■ ■ ■ MS78 




lf[ t ffl l IllTffTpT[ l 




mm 



T =5.0x10 K 



! \ 

9 * 



RM 

M85 

■ ■ ■ MS78 



<^ I liiirfWMHTfir 




ii^iii^n^rlT I liilltllMMtfiT I I I Iiiiiin>»iiijiiiyiii^ 



T =1.0x10 K 



DW 

RM 

M85 

■ ■ ■ MS78 



w 



/A 

iimttlLnt I I I I I I lii|||||i|Mii"^"tf 



7710 7725 7740 7755 7770 7785 7800 7815 7830 
Photon Energy (eV) 

Figure 2. Emission spectra of a low density {Ne = lO^^'cm"'') Ni plasma at three different temperatures in the high-temperature range 
where the G ratio from the models is in poor agreement. Spectra from the DW, M85, and MS78 models have all been red shifted by 9 
eV to facilitate an easier comparison. 



Moving on to a comparison of the effective recombina- 
tion rate coefficients in the four models (again presented only 
for Fe since the Ca and Ni results show similar trends) shown 
in figure [6] there are obvious discrepancies at high temper- 
atures and the cause of the differences in simulated spectra 
becomes evident. The two, more approximate models (MS78 
and M85) are in reasonable agreement with the more ac- 
curate DW and RM models in the middle-temperature re- 
gion. Outside the range, there is a marked difference (as 
much as a factor of five in the high-temperature region). 
This difference in the recombination rate coefficients in the 
high-temperature regime is responsible for the correspond- 
ing difference in the G ratio among the various calculations. 



The values of the new DW and RM recombination data are 
significantly smaller than the MS78 data, which result in de- 
creased recombination rates for populating the 'x', 'y', and 
'z' lines, while the impact-excitation rate for feeding the 'w' 
line remains essentially fixed for each type of calculation. 
Thus, the reduced DW and RM recombination data cause a 
delay in the onset of any increase in the G ratio to a higher 
temperature, resulting in the large differences displayed in 
figure [1] 

It should be noted that iNahar et al.l (|200ll ') saw a 
similar difference in the recombination rates in the high- 
temperature region when they compared their level-specific 
unified recombination rates with the direct component of 
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Figure 3. Emission spectra of a low density {Ne = lO^'^em^'') Fe plasma at three different temperatures in the high-temperature range 
where the G ratio from the models is in poor agreement. Spectra from the DW, M85, and MS78 model have all been red shifted by 7 eV 
to facilitate an easier comparison. 



the rec ombination rate expressions of iMewe &i Schriiveil 
(|l978al ). Again there is exceUent agreement between the DW 
and RM models over the entire temperature range, bolster- 
ing confidence in their validity. 



5 CONCLUSION 

While the results of calculations using recombination rates 
derived from more accurate and extensive distorted-wave 
and R-Matrix calculations agree quite well, both disagree 
significa ntly with the previous, mo r e approximate calcula- 
tions o f [Mewe fc Schriived (Il978a ): iMewe fc Gronenschijdl 
lIlQSll 'li lBelv-Dubau et all (|l982al lbl) and lMewe et al.l(|l985h . 



The results obtained from models based on the more accu- 
rate distorted-wave and R-Matrix data show a continually 
decreasing G ratio over the investigated temperature range. 
This behaviour disagrees strongly with the sharp upturn ob- 
served at high temperatures in the G ratios calculated from 
the more approximate models. The discrepancy is up to a 
factor of six in the G ratio (see figure[T]), and was determined 
to be caused by large differences in the recombination con- 
tributions. The more approximate rates produce results that 
agree reasonably well with the newer DW and RM data up 
to approximately 10* K, but can lead to unphysical values 
when used outside of this range. 

Further work is required to refine and verify calcula- 
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Figure 4. Emission spectra of a low density {N^ = lO^^'cm^'^) Ca plasma at three different temperatures in the high-temperature range 
where the G ratio from the models is in poor agreement. 



tions of the G ratio and the underlying data from which it 
is determined. Studies similar to the present one are needed, 
but including electron-impact excitation rates as well as col- 
lisional ionisation rates from the ground state of the Li-like 
species into the excited states of the He-like species. Inves- 
tigation of the effect of the ionisation balance may also be 
required. These studies are underway. 
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